Approximate black holes from a variable cosmological constant 
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The small or zero cosmological constant, A, probably results from a macroscopic cancellation mech- 
anism of the zero-point energies. However, nearby horizon surfaces any macroscopic mechanism 
is expected to result in imperfect cancellations. A phenomenological description is given for the 
residual variable cosmological constant. In the static, spherically symmetric case it produces ap- 
proximate black holes. The model describes the case of exponential decay by □ In A = —3a, were a 
is a positive constant. 

PACS numbers: 04.70.-s 
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The modern saga of the cosmological constant prob- 
lem concerns the divergence of the zero-point energies 
0. A particularly compelling mechanism for cancel- 
lation of the resulting cosmological constant has been 
proposed by Coleman in the Euclidean path integral ap- 
proach Q. Coleman's analysis predicts a zero expec- 
tation value of the cosmological constant by integration 
over thin but extended wormholes. However, a horizon 
surface is expected to perturb any such macroscopic can- 
cellation mechanism, by suppression of those wormholes 
which would tend to penetrate the surface. This pertur- 
bation suggests the possibility for a finite cosmological 
constant to reappear in the vicinity of horizon surfaces. 
In this Letter, a simple phenomenological approach is de- 
scribed to illustrate this scenario, using a self-consistent 
classical description of a variable cosmological constant. 

A space-time with finite cosmological constant, A, 
filled with a fluid with stress-energy tensor, T Q p, gives 
rise to the equations of motion 



G 



a/3 



'Ag a p + S a p = T a p. 



(1) 



We shall look for time-independent, spherically symmet- 
ric solutions in the Schwarzschild line-element 



As 2 



3 A dr 2 



r 2 (d6 2 +sin 2 



'At 1 . 



(2) 



A spherical region with finite A must be suitably closed 
by a shell with positive hoop stress to sustain the negative 
(radial) pressure inside. The shell also resides in a deep 
gravitational potential well, and therefore is expected to 
accumulate radiation and particles. If the particles inter- 
act infrequently (relative to the orbital time-scale), the 
stress tensor will develop an anisotropic pressure dom- 
inated by random orbital motion. In the limit of zero 
radial pressure of the fluid, its stress-energy tensor be- 
comes 

S a p = (E + P)u a Uj3 + Pk a/3 , (3) 

where vP = (0, 0, 0, e' v l 2 ) is the the average velocity 
four-vector of the particles, 



k af 3 = g a p - (d' r Ad 1 A)- 1 d a Ad A 



(4) 



is a projection tensor, E is the energy density of the fluid 
and P is its two-dimensional pressure. Hence, we have 



Tf = -diag(A,M, M,A + E), 



(5) 



were we have set M = A — P. It may be seen that 
the quantities A, E and P are scalars. The equations of 
energy-momentum conservation, 
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d a Tf 



■pa rpj 

1 ar/ ± p 
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reduce to 



(6) 



(7) 



With the tabulated expressions for the connections (e.g.: 
||), this obtains 



A' = -(M-A) + ^-E. 
r 2 



(8) 



At this point, recall the the weak energy condition on the 
stress-energy tensor 



A + E>{), 



(9) 



were the infimum is taken over all time- like vectors . 
Clearly, (||) allows for the fluid to consist not only of 
particles (E p > 0), but also of negative mass holes (Eh < 
0) in regions with A > 0. Thus, coexistence of particles 
and holes makes it possible to have positive hoop stress, 
even with E = 0. 

We will look for Schwarzschild solutions which join to 
space-time with a variable cosmological constant in the 
neighborhood of the origin, that is 



1 - ^- , e u ~ 1 - E, P ~ 0, (10) 



were Ao is the value of the cosmological constant at the 
origin. The relevant Einstein equations are those for G r r 
and G\, 



1 



e~ A (^ + ^)-^ = -A, 
-XfX' 
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-(A + E). 



These may be written as 



rv' = -1 + e A (l - r 2 A), 
rX' = l-e x (l-r 2 [A + E}). 



(11) 



(12) 



Whenever A = E = P = beyond some r = tq, ( |l2| ) 
recovers the Schwarzschild solution, though perhaps after 
a rescaling of the time coordinate (in case of X + v ^ 0) to 
maintain the familiar relationship —gttgrr — 1- Clearly, 
the simplest case follows with E = 0, which is allowed 
in our particle- hole model if any gravitational separation 
of particles and holes can be neglected on the scale of 
variation in A. 

The simplest case, therefore, is given by a constitutive 
relation only for the anisotropy, P, in the pressure. The 
regularity condition of P at the origin, P = 0(r 2 ) for 
small r, translates covariantly to 

(13) 

F ~ F(X) ' [16 > 

were F (A) is a continuous function of A. For example, 
F(X) = 2aA a , so that 

'~^(A') 2 (14) 



P = 



2aA 

giving a reduction of (||) to 

A' = ~are x A a 



(15) 



Equation (fTq) shows a positioning of a transition layer in 
A at large er, i.e., at a surface of large red shift. In the 
case of a = 1, the covariant generalization of (Il5|) is 



□ In A = -3a. 



(16) 



A solution is illustrated in Figure 1. With < a < 1, the 
distinguishing feature arises of A vanishing at finite r, as 
illustrated in Figure 2. The covariant generalization of 
(|l|) for a ^ 1 is given by 

□A (1 ~ Q) = -3a(l-a); (17) 

for a > 1, A(r) decays algebraically in the far field limit. 



F(A)=2aA 




Figure 1. Distributions of the metric quantities A, v and 
A as a function of radius in the case of F = 2aA (a = 1) for 
a = 0.5, 1 and A = 5 at the origin. A high red shift accom- 
panies the transition layer in A. Outside of this layer, the 
cosmological constant is exponentially small, and space-time 
is essentially described by the exterior Schwarzschild solution. 
Inside this layer, the cosmological constant is effectively con- 
stant. The origin is singulaty free. 



F(A)=2aA u 
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Figure 2. Distributions of the metric quantities A, v and A 
as a function of radius in the case of F = 2aA ' 5 (a = 0.5) 
for a = 1, 2 and A = 5 at the origin. The transition layer is 
sharper compared with the previous case F = 2aA with the 
same a = 1 as displayed in Figure 1. Furthermore, a bifur- 
cation surface appears beyond which A vanishes identically, 
slightly outside the surface of high red shift. The origin is 
singularity free. 

It is an open question whether ( |l7|) preserves non- 
negativity of A' 1- ") (Ao > 0) in dynamical space-times. 
It may take numerical simulations to study the behavior 
of the surface of bifurcation (at which A becomes zero) 
in this setting. For this reason, a preliminary preference 
may be given to (|l^) which has a — 1. It should be 
metioned that any terms on the right hand-side of ( [Trjf - 
|l7| ) which vanish for E = have been suppressed. 

The existence of solutions can be seen from bounded- 
ness of A. (|l2|) shows that A can blow up at a finite radius 
only when r^A > 1. By @, rX' > -(e A - 1), so that 
A > 0. Then A' < re A A by the same equation. By ([IF 
re A A can be expressed in terms of A', so that 



A' < - 



1 A' 



A«- 



It follows that 



w x w x 1 /" A(r0) dA 

A(r) < A(r + - 
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Clearly, if < a < 2, we may take r = 0, so that 



(18) 



(19) 



2 



< A < 



a(2-a) 



AI — ol , 
n < OO. 



(20) 



For a > 2, choose ro such that A(ro)rQ > 1. For those 
T > ^0; where A(r)r 2 > 1, we have 



A(r) < A(r ) + i/ 1 A / ( ; 2 o) A 1 -^A 



< A(r ) + i 



lnA(r )r 2 (a = 2), 
r 2 ( Q ~ 2 ) (a > 2). 



(21) 



Again, no blow-up at finite radius occurs. 
The fully covariant stress-energy tensor is 



S a p = Eu a Uf3 + Ph a/ 3 



dgAdpA 



(22) 



with h a p = g a p + u a up. The energy E has been included 
for consistency in view of the normalization u 2 = — 1 , so 
that the four equations of energy-momentum conserva- 
tion, 



V a 7> Q = 



(23) 



can be satisfied with the three degrees of freedom in u 13 
plus the one degree of freedom in E. As we have seen, E 
will depart from zero only in genuinely dynamical space- 
times. The evolution equation for E and u@ satisfies 



V a {Eu a up) + V a S% - dp A = 0; 



in particular, 



V a {Eu a ) = u^dpA 



(24) 



(25) 



expresses the evolution of E. 

The present solutions feature surfaces of arbitrarily 
high red shift, outside of which space-time is essentially 
that of general relativity. These solutions are related to 
the earlier solution of van Putten M, which were pro- 
posed as as approximate black holes for numerical reg- 
ularization. Spherically symmetric static solutions with 
surfaces of arbitrarily high red shift have recently also 
been described by 't Hooft ||||, in an entirely different 
context. These solutions are both hampered by naked 
singularities at the origin (though of different kinds). Nu- 
merical simulations by Coptuik et al. Q also show the 
solution of to be unstable. This naturally raises the 
same question of stability of the present solutions, some- 
thing which falls outside the scope of the present work. 
It should be emphasized that the present solutions are 
an improvement (over foremcntioncd related solutions) 
in that these are free of any singularities. 

In regards to the application to numerical relativ- 
ity, note that the stress-energy tensor remains differen- 
tiable everywhere, though with a discontinuity in its first 
derivative if < a < 1. This still provides a proper 
source term in the covariant, strictly hyperbolic formula- 
tion in terms of nonlinear wave equations by van Putten 
& Eardley , given by 



Ouj aab - W a uj lab - [w 7 , V a uj 7 } ab = 167TT Qah , (26) 

were u> aa b are the tetrad connections of the tetrad ele- 
ments (e a ) a (here, Latin indices refer to the tetrad in- 
dices and Greek indices refer to the coordinates). The 
source term T aab = r ct/37 (e a )' 9 (e( ) ) 7 , is given by 



(27) 



were T = T 7 is the trace of the stress-energy tensor. The 
equations in the exterior vacuum, therefore, are 



OuJ a ab ~ [^ 7 , V Q U/ 7 ] a6 = 0. 



(28) 



These have been given a second-order numerical imple- 
mentation in one dimension || . In the case of < a < 1 , 
T aa b in ( p6| ) is at most discontinuous, due to second 
derivatives of A; for a > 1, there are no discontinuities. 
It is expected, therefore, that for all a > the stress- 
energy tensor (^), plus —Ag a p and constitutive relation 
( |l3| ) , can indeed be integrated into this approach of nu- 
merical relativity. 
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